The universality class of absorbing phase transitions with a conserved field 
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We investigate the critical behavior of systems exhibiting a continuous absorbing phase transition 
in the presence of a conserved field coupled to the order parameter. The results obtained point out 
the existence of a new universality class of nonequilibrium phase transitions that characterizes a 
vast set of systems including conserved threshold transfer processes and stochastic sandpile models. 
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Absorbing phase transitions (APT) are a category of 
critical nonequilibrium phase transitions, widespread in 
condensed matter physics and population and epidemics 
modeling ffl. Directed percolation (DP) [QJ|] has been 
recognized as the paradigmatic example of a system ex- 
hibiting a transition from an active to a unique absorbing 
phase. DP defines a precise universality class (theoreti- 
cally described by the Reggeon field theory |],|) which 
has proven to be very robust with respect to the intro- 
duction of microscopic modifications. The Reggeon field 
theory is at the heart of a strong claim of universality, 
summarized in the following conjecture [pi: Continuous 
absorbing phase transitions to a unique absorbing state 
fall generically in the universality class of directed per- 
colation. This conjecture is expected to hold for models 
with short range interactions that, most importantly, do 
not posses additional symmetries. 

Many examples of APT subject to extra symmetries, 
and thus out of the DP class, have been identified in re- 
cent years. Among them we find systems with symmet- 
ric absorbing states || , models of epidemics with perfect 
immunization (the so-called dynamic percolation class) 
0, and systems with an infinite number of absorbing 
states H . Very recently, it has been pointed out that the 
critical point of self-organized critical (SOC) sand- 
pile models can also be interpreted as a continuous phase 
transition with many absorbing states [[ll|]l2[ | . What dis- 
tinguishes sandpile models from other models with ab- 
sorbing states, is that the control parameter, represented 
by the global density of particles, is a conserved quantity. 

Given the large class of systems whose dynamics in- 
volves conserved fields, it becomes particularly interest- 
ing to explore in general the effect of conservation rules in 
APT. With this purpose in mind, in this Letter we report 
the critical behavior of several models showing absorbing 
transitions that strictly conserve the number of particles 
or energy. In particular, we introduce a conserved lattice 
gas (CLG) [[13| with short range stochastic microscopic 
dynamical rules, that undergoes a continuous phase tran- 
sition to an absorbing state at a critical value of the par- 
ticle density. We present extensive numerical simulations 
in d = 2 of the stationary and spreading properties of the 
model, and determine the full set of critical exponents. In 
order to prove definitively the existence of a well-defined 
universality class we have also performed simulations of 



a conserved threshold transfer process (CTTP) jTJ], and 
several fixed energy sandpile models with stochastic rules 
[fTo] |l2| . All models provide critical exponents compat- 
ible with a single and broad universality class that em- 
braces all APT in stochastic models with a conserved 
field. This evidence leads us to conjecture that, in ab- 
sence of additional symmetries, absorbing phase transi- 
tions in stochastic models with infinite absorbing states 
and activity coupled to a static conserved field define a 
unique and per se universality class Jl5| . This result is 
relevant in the understanding of several reaction-diffusion 
systems, sandpile models and activated processes that 
could share the same theoretical description. 

The CLG model is defined on a d-dimcnsional square 
lattice. To each site i it is associated a binary variable 
ni that assumes the values rn = 1 if the site is occu- 
pied by a particle or n, = if the site is empty. Double 
occupancy is strictly forbidden. Nearest neighbors par- 
ticles repel each other via repulsive short range interac- 
tions. As a product of this interaction, at each time step 
particles with nearest neighbors jump into one of their 
empty nearest neighbor sites, selected at random. The 
only dynamics in the model is due to these active parti- 
cles; isolated particles do not move. The dynamics can be 
implemented with either sequential or parallel updating. 
In the latter case, an exclusion principle is applied so that 
two particles never attempt to move into the same site. 
We impose periodic boundary conditions, and since the 
dynamics admits neither input nor loss, the total num- 
ber of particles N — J2i n i(t) * s a conserved quantity. 
It is clear that the model allows an infinite number (in 
the thermodynamic limit) of absorbing configurations, in 
which there are no nearest neighbor particles. 

In the CLG model, the constant particle density n = 
N/L d acts as a tuning parameter. Initial conditions are 
generated by placing at random in the lattice nL d parti- 
cles, generating an homogeneous and uncorrelated distri- 
bution. For small densities, the system will very likely fall 
into an absorbing configurations with only isolated parti- 
cles. For large densities, the system reaches a stationary 
active state with everlasting activity (this is trivially the 
case for n > 1/2). We shall see in the following that as 
we vary n, the CLG model exhibits a continuous transi- 
tion separating an absorbing phase from an active phase. 
The phase transition occurs for a nontrivial density n c 
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FIG. 1. Active-site density in surviving trials at the crit- 
ical point for the CLG. From top to bottom, system sizes 
are L = 160, 192, 256, 320, 416, and 512. Inset: Stationarity 
active-site density as a function of L for a critical (center), 
subcritical (down), and supercritial (up) particle density. 
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FIG. 2. Collapse of the stationary density p a (An, L)L' 3// " J - 
as a function of AnL 1 ^" 1 - for the CLG. The slope of the line 
is 0.63. 



(< 1/2). APT are characterized by the order parameter 
p a measuring the density of dynamical entities, in our 
case the density of nearest neighbor particles. The or- 
der parameter is null for n < n c , and follows a power 
law p a ~ (n — nc) 13 for n > n c . The system correlation 
length £ and time r both diverge as n —> nj. In the 
critical region the system is characterized by power law 
behavior, namely £ ~ (n — n c )~ u± - and r ~ (n — n c ) _,y ii . 
The dynamical critical exponent is defined as r ~ £ 2 , 
with z = f| | / . These exponents fully define the critical 
behavior of the stationary state of the model. 

In order to study the critical point of the CLG model, 
we performed numerical simulations in d = 2 for systems 
with size ranging from L — 64 to L — 512, averaging over 
10 4 — 10 5 independent initial configurations. Very close 
to the critical point we have £ 3> L, so that the actual 
characteristic length of the system is the lattice size L. 
Because of its finite size, the system will enter sometimes 
an absorbing configuration even for values of n in the 
supercritical region. It is then convenient to introduce 
averages over a set of independent trials and calculate 



the quasi-stationary properties in the active phase from 
a restricted average over surviving trials with nonzero 
final activity. 

As shown in Fig. ^, after a transient which depends 
on the system size L and An = n — n c , the surviving 
samples average of the density of active sites reaches a 
stationary state p a (L, An). Close to the critical point, 
the finite size scaling ansatz tells us that all quantities 
depend on the system size through the ratio L/£, and 
the order parameter follows the finite size scaling form 

n 

p a (An,L) = L-P/^g(LV^An) , (1) 

where Q is a scaling function with Q(x) ~ x@ for large x. 
For An = the stationary density follows the pure power 
law behavior p a ~ L~P/ U± . On the other hand, for values 
of n in the supercritical regime p a should be independent 
of L for L while in the subcritical regime p a should 
decay faster than a power law. This allows us to locate 
the critical value n c of the particle density as the only 
value of n at which we recover a nontrivial power law scal- 
ing for the density of active sites. In Fig. [I] we observe 
power law scaling for n = 0.23875, but clearly not for 
0.2387 or 0.2388, indicating that n c = 0.23875(5) (Fig- 
ures in parenthesis indicate the statistical uncertainty in 
the last digit). From the power law decay we find the 
exponent ratio f3/v± — 0.81(3). An independent esti- 
mate of the exponent f3 can be obtained by looking at 
the scaling of the active-site density with respect to An 
for the size L — 320. The resulting power law behav- 
ior yields (3 = 0.63(1), where the error is mainly due 
to the uncertainty in the critical point n c . A consis- 
tency test can be performed by considering the active 
site density away from the critical point. In Fig. || we 
plot p Q (An,L)L /5 /^ versus AnL 1 ^ for vy_ = 0.78, 
/3/V_i_ = 0.81 and n c = 0.23875. As one would expect 
all the data collapse onto a single curve, following the 
scaling form Eq. (lj). A further check is provided by the 
direct fitting of the large x behavior of the scaling func- 
tion Q{x) that gives [3 — 0.63, recovering the independent 
measurement at L — 320. 

To determine the dynamical exponents we turn our at- 
tention to the scaling properties of time dependent quan- 
tities. In particular, we can define a characteristic time 
by studying the decay of the probability Pit) that a ran- 
dom initial configuration has survived up to time t. At 
the critical point this probability decays, at large times, 
as P(t) ~ exp(— t/r). At An = the effective charac- 
teristic length is the system size L, and we have that 
t(L) ~ L z . We can access the value of t(L) by a direct 
fitting of the P(t) exponential tail; z is then estimated 
from the behavior of t(L) for different L. Again, a clean 
power law behavior is obtained for n c = 0.23875, yielding 
z = 1.52(6). Also in this case a consistency check can be 
performed by studying the time decay of the active sites 
density p a .aii(t), averaging over all trials, even those that 
have reached an absorbing state. Assuming a single char- 
acteristic time scaling as L z , we have at An = |ll| 
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Pa,all{t,L)=t- e T{tL- Z ) , (2) 

where F(x) is a constant for i < 1, and decays faster 
than any power law for x 3> 1. Data from simulations 
with different L can be collapsed onto a universal curve 
by plotting p a , a ii{t, L)t B versus tL~ z . The best collapse is 
obtained for 9 = 0.43 and z — 1.52, confirming the value 
obtained for the dynamical critical exponent. The expo- 
nent 9 = 0.43(1) is recovered also from a direct fitting of 
the decay of the stationary density averaged over surviv- 
ing trials (see Fig. |l|). In usual APT, the latter exponent 
obeys 9 = fi/v\\- This relation assumes a standard scal- 
ing behavior at An = for p a (t). In our model, however, 
the simple scaling behavior is broken by an anomalous 
scaling regime (visible in Fig. [l] by the sharp drop just 
before the stationary state) that seems to grow steeper 
with increasing L. It follows that data collapse in time 
is not achievable with standard scaling forms, and that 
9 violates the usual scaling relation. Albeit its origin is 
not clear, it is noteworthy that this anomaly is common 
to all APT with conserved fields inspected so far ]n| , |l2"[ . 
At this respect, it is interesting to note that the expo- 
nents of the model fulfill all scaling relation in standard 
APT. In what respects to hyperscaling relations, some of 
them, as D = d + z — P/v±, are fulfilled, while others, 
like rj + 5 + 9 = d/z, arc not. This is again due to the 9 
exponent anomaly. 

In APT it is possible to obtain more information on 
the critical properties by studying the evolution (spread) 
of activity in systems which start close to an absorbing 
configuration []3j. In each spreading simulation, a small 
perturbation is added to an absorbing configuration. It is 
then possible to measure the spatially integrated activity 
N(t), averaged over all runs, and the survival probability 
P(t) of the activity after t time steps. Only at the critical 
point we have power law behavior for these magnitudes. 

Here, we will follow a procedure equivalent to the def- 
inition of slowly driven simulations in sandpiles, that en- 
lights the connections with these models. Instead of fix- 
ing the density n by working with periodic boundary con- 
ditions, and thus studying the system at a given distance 
below the critical point, we impose open boundary con- 
ditions and start each spreading experiment by adding a 
new particle. Under these conditions, the system flows to 
a stationary state with balance between the input of par- 
ticles and the boundary dissipation. In the limit in which 
the particle addition is infinitely slow with respect to the 
spreading of activity, the system reaches a critical state 
with density n c (in the thermodynamic limit) fl7|| . The 
infinitely slow drive is implemented by adding a new ac- 
tive particle only when the system falls into an absorbing 
configuration. The system thus jumps between absorb- 
ing states via avalanche-like rearrangements, and we can 
associate each spreading experiment with an avalanche. 
The probability distribution P s (s) of having a spread- 
ing event involving s sites, as well as the the quantities 
N(t) and P(t) can be measured. The only characteristic 
length is the system size L, and we can write the scaling 







Steady state exponents 
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CLG 
CI IP 

Manna 
DP 


0.63(1) 
0.64(1) 
0.64(1) 
0.583(4) 


0.81(3) 1.52(6) 0.43(1) 
0.78(3) 1.55(5) 0.43(1) 
0.78(2) 1.57(4) 0.42(1) 
0.80(1) 1.766(2) 0.451(1) 








Spreading exponents 
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CLG 
CTTP 
Manna 
DP 


1.29(1) 
1.28(1) 
1.28(1) 
1.268(1) 


2.75(1) 1.53(2) 0.29(1) 
2.76(1) 1.54(1) 0.30(3) 
2.76(1) 1.55(1) 0.30(3) 
2.968(1) 1.766(2) 0.230(1) 


0.49(1) 
0.49(1) 
0.48(2) 
0.451(1) 



TABLE I. Critical exponents for spreading and steady 
state experiments. Figures in parenthesis indicate the sta- 
tistical uncertainty in the last digit. Steady state Manna ex- 
ponents from Ref. 



forms H: 

N(t) = f'f(t/L z ), P(t) = t- s g(t/L z ), 
P s ( S ) = S ~ T °h( S /L D ), (3) 

where the scaling functions f(x), g(x) and h(x) are de- 
creasing exponentially for x S> 1, and we have considered 
that the spreading characteristic time and size are scaling 
as L z and L D , respectively. Simulations were performed 
for system of size between L = 64 and L = 1024, averag- 
ing over at least 5 x 10 6 spreading experiments. The ex- 
trapolation of the measured densities at infinite L yields 
a critical density n c = 0.2388(1), in perfect agreement 
with steady state simulations. The scaling exponents 
are measured using the now standard moment analysis 
technique ]lqjl9| . The resulting exponents are summa- 
rized in Table I. In particular, the dynamical exponent 
z = 1.53(2) is in excellent agreement with the station- 
ary state simulations, confirming the presence of a single 
critical behavior for both cases. 

In order to provide further evidence for the existence 
of a general universality class, we have simulated several 
other models exhibiting an APT in the presence of a con- 
served field. The first is a conserved threshold transfer 
process (CTTP). In the CTTP, the sites of a lattice can 
be vacant, singly occupied, or doubly occupied by parti- 
cles, corresponding to a dynamic variable rij = 0, 1 or 2, 
respectively. Values n, > 2 are strictly forbidden. Dy- 
namics affects only doubly occupied sites: every site with 
rii = 2 tries to transfer both its particles to randomly se- 
lected nearest neighbors with rij < 2. Singly occupied 
sites are, on the other hand, inert. The total number of 
particles N — rii is thus constant in time. Results 
from simulations are obtained along the lines shown for 
the CLG and are reported in Table |; in this case, the 
largest sizes used are L = 512 for the stationary expo- 
nents and L = 1024 for the spreading exponents. As an 
example of our simulations, in Fig. |^ we plot p a _ a u (£, L)t e 
as a function of tL~ z , which shows a remarkable data 
collapse. We have also investigated the Manna sandpile 
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FIG. 3. Collapse of the active-site density p a ,aii(t, L)t as 
a function of tL~ z for the CTTP. System sizes range from 
L = 64 to L = 384. 



model, and its variations by the inclusion of a stochastic 
threshold pC(| . In this case, an absorbing phase transi- 
tion is obtained by using periodic boundary conditions 
and a fixed number of sand grains (energy) as reported 
in (ll],[l2| ■ All models lead invariably to the same univer- 
sality class as the Manna model |2]J] (complete results on 
these models will be reported elsewhere). 

Our results provide striking evidence for a common 
critical behavior which is incompatible with the DP uni- 
versality class (see Table [j]) plfl . More noticeably, the 
models share also the same scaling anomaly in the ex- 
ponent 6 ]l2| ], signalling a common behavior in the tran- 
sient regime to the stationary state in the active phase. 
This uniformity of results confirms the hypothesis of a 
unique universality class for all models with the same 
conservation symmetry, and lead us to conjecture that, 
in the absence of additional symmetries, absorbing phase 
transitions in systems with stochastic dynamics in which 
the order parameter is locally coupled to a static con- 
served field define a single and new universality class. 
This conjecture is supported by noticing that the sandpilc 
model has the same structure and basic symmetries of the 
present CLG, once the field n(x, t) is replaced by local 
energy (sand grains) field Indeed, in all models 

presented here, a conserved non-critical field is dynami- 
cally coupled to the non-conserved order parameter field 
p a - Very likely this basic structure will be reflected in 
a unique theoretical description (a field theory with the 
same relevant terms and symmetries) that accounts for 
the shared critical properties of these models. 
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